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Abstract 


Background: Infectious bronchitis virus (IBV) is first to be discovered coronavirus which is probably endemic in all 
regions with intensive impact on poultry production. In this study, we used two-dimensional gel electrophoresis (2- 
DE) and two-dimensional fluorescence difference gel electrophoresis (2-DIGE), coupled with matrix-assisted laser 


desorption/ionization time-of-flight tandem mass spectrometry (MALDI-TOF/TOF-MS), to explore the globa 
es of trachea and kidney tissues from chicken at different stages infected in vivo with the highly 


proteome profi 


virulent ck/CH/LDL/97I Ps strain of infectious bronchitis virus (IBV) and the embryo-passaged, attenuated c 


LDL/97I P15 strain. 


CH/ 


Results: Fifty-eight differentially expressed proteins were identified. Results demonstrated that some proteins which 
had functions in cytoskeleton organization, anti-oxidative stress, and stress response, showed different change 
patterns in abundance from chicken infected with the highly virulent ck/CH/LDL/97I Ps strain and those given the 
embryo-passaged, attenuated P1;5 stain. In addition, the dynamic transcriptional alterations of 12 selected proteins 
were analyzed by the real-time RT-PCR, and western blot analysis confirmed the change in abundance of heat 


shock proteins (HSP) beta-1, annexin A2, and annexin AS. 


Conclusions: The proteomic alterations described here may suggest that these changes to protein expression 
correlate with IBV virus’ virulence in chicken, hence provides valuable insights into the interactions of IBV with its 
host and may also assist with investigations of the pathogenesis of IBV and other coronavirus infections. 


Keywords: Infectious bronchitis virus, Proteomics, Chicken, Trachea, Kidney 


Background 

Coronaviruses (CoVs) are enveloped single-stranded 
positive sense RNA viruses that belong to the family 
Coronaviridae in the order Nidovirales. They are able to 
infect humans as well as other animals, including cows, 
pigs, mice, and chickens, they generally cause respiratory 
infection, gastrointestinal, and neurological disorders of 
varying severity. Infectious bronchitis virus (IBV) was 
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the first coronavirus to be discovered, and is classed 
among the Gamma coronaviruses on the basis of anti- 
genic and genetic relatedness [1]. It is a major poultry 
pathogen and is probably endemic in all chicken-raising 
regions; it has a severe impact on poultry production, 
causing heavy economic losses. All strains of IBV are 
capable of infecting a large range of epithelial surfaces 
of chickens, such as those of the trachea, kidney, oviduct 
and proventriculus [2]. 

Coronavirus infection has dramatic effects on host cell 
morphology, transcription and translation patterns, the 
cell cycle, cytoskeleton, suppression of interferon, and 
apoptosis pathways. Coronavirus infection may also 
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cause inflammation, alter the immune and stress 
responses, and modify the coagulation pathways [3]. 
Such profound functional and morphological changes in 
host cells are associated with significant changes in the 
patterns of expression of host cell genes. Several studies 
have described changes in host gene expression asso- 
ciated with coronavirus infection, as documented by 
microarray technologies [4-9]. But ultimately, protein 
expression and post-translational modification (PTM) 
determine virus replication. Furthermore, transcriptome 
analyses only provide a snapshot of gene expression pat- 
terns; they also suffer from several limitations, including 
inconsistencies with the levels of expression of the cor- 
responding proteins as well as lacking the ability to pro- 
vide information on PTM. Approaches that use 
proteomics are promising because they can circumvent 
some of the issues associated with transcriptomics 
approaches [10]. More recently, comparative proteomics 
analysis has emerged as a valuable tool for the establish- 
ment of the global host protein profile in response to 
virus infection. It has been used to study enveloped 
RNA viruses such as influenza virus, respiratory syncy- 
tial virus (RSV), parainfluenza virus (PIV), human 
metapneumovirus (hMPV), SARS-CoV, and mouse 
hepatitis virus (MHV) [11-18]. It provides invaluable 
information on the cellular signaling pathways involved 
in either the cellular response to viral infections, or the 
viral manipulation of cellular machinery to ensure their 
own survival. For IBV, to the best of our knowledge, 
only some recent studies have investigated the changes 
in the expression of cellular proteins during IBV infec- 
tion in ex vivo or in ovo [19-21]. However, the in vivo 
infection model could yield more biologically relevant 
insights into pathogenesis. 

In this study, we used two-dimensional gel electro- 
phoresis (2-DE) and two-dimensional fluorescence dif- 
ference gel electrophoresis (2-DIGE), coupled with 
matrix-assisted laser desorption/ionization time-of-flight 
tandem mass spectrometry (MALDI-TOF/TOF-MS), to 
explore global changed proteome profiles of trachea and 
kidney tissues from chicken at different stages infected 
in vivo with the highly virulent ck/CH/LDL/97I Ps strain 
of IBV and an embryo-passaged strain of attenuated 
virulence, ck/CH/LDL/97I P 45. In total, 58 differentially 
expressed proteins were identified and classified into 
several functional categories, including cytoskeleton 
organization, anti-oxidative stress, the stress response, 
acute phase response, and energy metabolism. In addi- 
tion, the dynamic transcriptional alterations of 12 
selected proteins were analyzed by the real-time RT- 
PCR method. Simultaneously, western blot analysis con- 
firmed the change in abundance of the heat shock pro- 
teins (HSP) beta-1, annexin A2, and annexin A5. The 
potential roles of some of these identified proteins are 
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discussed in order to characterize their potential func- 
tional roles during IBV infection in vivo. These results 
provide valuable insights into the interactions of IBV 
with its host, and may also be useful in investigations of 
the pathogenesis of IBV and other coronaviruses. 


Results 
IBV antibody detection and observed clinical signs 
All chickens exhibited respiratory clinical signs at about 
4-14 dpi with the IBV ck/CH/LDL/97I Ps strain. The 
clinical signs included tracheal rales, watery eyes, nasal 
mucus, and sneezing. The clinical signs shown by the 
inoculated birds tended to disappear gradually after 14 
dpi. Gross lesions of the chickens killed in the P;- 
infected group were confined mainly to the kidneys. The 
kidney parenchyma of the dead birds was pale, swollen 
and mottled; the tubules and urethras were distended 
with uric acid crystals [22]. For chickens in the control 
group and those inoculated with the IBV ck/CH/LDL/ 
971 Pi15 strain, no respiratory clinical signs and no gross 
lesions were observed during the experimental period. 
As summarized in Additional file 1, no chickens 
inoculated with the ck/CH/LDL/97I Ps; or the ck/CH/ 
LDL/97I P,15 strain of IBV showed seroconversion at 4 
dpi. Antibodies appeared at 7 dpi and all of the chickens 
showed seroconversion after 14 dpi with both IBV 
strains. The chickens in the non-inoculated control 
group showed a negative serum antibody response. 


Analysis of viral load in the trachea and kidney of IBV- 
infected chickens 

Successful IBV infection was also verified using real- 
time RT-PCR. Results are presented in Figure 1. Virus 
was not detected from the trachea and kidney of chick- 
ens from control group. In the Ps-infected group, virus 
was detected from trachea and kidney at 4, 7, 14, and 
21 dpi, the peak of viral copy number was reached at 4 
dpi, after which time viral load fell. In trachea of P,5- 
infected group, the peak of viral copy number was at 4 
dpi, then fell until 14 dpi, and a little fluctuated at 21 
dpi. In kidney of P,,5-infected group, the peak of viral 
copy number was also at 4 dpi, then fell at 7 dpi, and 
only a little fluctuated at 14 and 21 dpi. Furthermore, 
both in trachea and kidney, the P,,5-infected group had 
the lower viral genome copies than Ps-infected group. 


Analysis of differentially changed proteins in abundance 
by 2-DE and 2-DIGE 

The global protein changes in the chicken trachea and kid- 
ney tissues at different stages after infection with the IBV 
ck/CH/LDL/97I Ps and ck/CH/LDL/97I Pj;5 strain were 
investigated. For the tracheal protein samples, Figure 2 
shows representative images of the tracheal samples: 1366 
+ 39, 1536 + 126, and 1600 + 167 protein spots were 
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Figure 1 The viral load in each sample was quantified using real-time RT-PCR. The average viral copy number (VCN) per g tissue of each 
group was calculated. Error bars indicate standard error of the mean, and dpi represent days post-inoculation. 


detected in gels from the control group, P;-infected group, 
and P,,5-infected group at 4 dpi; 1355 + 300, 1518 + 175, 
and 1078 + 122 protein spots were detected in gels from 
the control group, Ps-infected group, and P,,5-infected 
group at 7 dpi; 1293 + 91, 1365 + 126, and 1220 + 56 pro- 
tein spots were detected in gels from the control group, 
P;-infected group, and P,,5-infected group at 14 dpi; 1204 
+ 91, 1236 + 42, and 1111 + 50 protein spots were 
detected in gels from the control group, P;-infected group, 
and P,,5-infected group at 21 dpi. For the kidney protein 
samples, Figure 3 shows representative images of the kid- 
ney samples at 4, 7, 14, and 21 dpi: 2315 + 87, 2482 + 189, 
2607 + 238, and 2593 + 192 protein spots were detected, 
respectively. 

The number of differentially expressed protein spots 
in the Ps and P,,5 infected groups is summarized in 


Additional file 2. In the tracheal total proteome at 4 dpi, 
the host response to infection with both Ps and Py15 
strains appears to involve predominantly increase of 
proteins abundance. By 7 dpi, the pattern is different, 
with infection with the Pj;5 strain resulting in a predo- 
minant increase of proteins abundance, and with P; 
exhibiting a majority of decrease events. By 14 and 21 
dpi, both viruses produce a broadly similar response, 
with the vast majority of changes in protein abundance 
involving increase. For the kidney total proteome, at 4 
dpi, Ps infection induced an approximately equivalent 
number of proteins with increased and decreased abun- 
dance, whereas P;,5 infection induced mainly increase 
of proteins abundance. By 7, 14, and 21 dpi, both P; 
and Py,5 infection induced mainly increase of host pro- 
teins abundance. 
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Figure 2 Representative 2-DE maps of tracheal tissues from the IBV ck/CH/LDL/97I Ps-infected group, P115-infected group and control 
group. (A) 4 dpi; (B) 7 dpi; (CQ) 14 dpi; (D) 21 dpi. Protein samples were separated on 13 cm pH 3-10 linear IPG strips, followed by SDS-PAGE, 
and stained with Coomassie Blue R-350. The images were analyzed using Image Master 2D Platinum 6.0 software. The differentially expressed 
protein spots identified were marked with circles and labeled with the respective Match ID listed in Tables 1 and 2. 


Identification of differentially changed proteins in TOF-TOF MS. From trachea and kidney tissues, 24 and 
abundance by MS 34 proteins were successfully identified, respectively. 
All interesting changed protein spots in trachea and kid- Detailed information on the identified proteins is pro- 
ney tissues were excised, and identified by MALDI- vided in Tables 1, 2, 3 and 4, Additional file 3, and 
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Figure 3 Representative 2-DIGE maps of kidney tissues from the IBV ck/CH/LDL/97I P;-infected group, P,,5-infected group and control 
group. (A) 4 dpi; (B) 7 dpi; (CQ) 14 dpi; (D) 21 dpi. Protein samples labeled with fluor were separated on 24 cm pH 3-10 linear IPG strips, followed 
by SDS-PAGE, gels were scanned on a Typhoon 9400 scanner, and image analysis was performed with Ettan™ DeCyder Software version v6.5. 


The protein spots identified were marked with circles and labeled with the respective Match ID listed in Tables 3 and 4. 
TY 


Table 1 Similar abundance changed proteins in tracheal from chicken infected with IBV ck/CH/LDL/971 P; and ck/CH/LDL/97I P1415 


Match Protein description Accession Protein Protein Mw pl Average ratio and p-value 
ID? no.” score score Cl (Da) 
%* 
P;-infected/control P115-infected/control 
4 dpi 7 dpi 14 dpi 21 dpi 4 dpi 7 dpi 14 dpi 21 dpi 
ratio p ratio p_ ratio p_ ratio p_ ratio p_ ratio p_ ratio p_ ratio p 
4467 annexin Al (ANXA1) gi 910 00 38761 7.05 0.89 0.435 049 0.001 
46195459 
2051 B-creatine kinase gil211235 631 00 42525.5 5.78 144 0.025 1.77 0.006 
926, cartilage matrix protein gi 314 00 54798.1 663 2.71 0.006 1.24 0427 194 0.003 248 0.025 1.96 0.012 246 0 
736, precursor 71896317 
817, 
4581 
987 fibrinogen beta chain gil211780 612 00 532718 7.18 2.75 0.006 2.57 0.011 
4679 glutathione S-transferase gil4959550 79 99.958 213265 8.03 136 0.005 18 O 
class-alpha 
4036 keratin, type | cytoskeletal gi 772 00 46111 4.94 243 0.016 2.15 0.009 
19 45384356 
5150 nmrA-like family domain- gi 653 00 32726 7.77 221 0.004 1.72 0.009 
containing protein 1 71897147 
4405 PREDICTED: microfibril- gi 331 00 28600 5.2 0.57 0.001 0.97 0.695 
associated glycoprotein 4. 18470900 
(EST) 
2021 Vimentin (VIM) gil212868 1160 00 531666 5.09 0.76 0.015 0.56 0.001 


a) Match ID represents unique numbers assigned to each spot in the matched standard of ImageMaster 2D Platinum version 6.0 software 
b) Accession no. is the MASCOT result of MALDI-TOF/TOF searched from the NCBInr database 
c) The criterion for a successfully identified protein is a protein score confidence interval (C.l. %) for PMF and MS/MS data = 95% 
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Table 2 Differential abundance changed proteins in tracheal from chicken infected with IBV ck/CH/LDL/97I P; and ck/CH/LDL/97I P,45 


Match Protein description Accession Protein Protein Mw p/ Average ratio and p-value 
ID? no.” score score Cl (Da) 
%S 
P.-infected/control P,,5-infected/control 
4 dpi 7 dpi 14 dpi 21 dpi 4 dpi 7 dpi 14 dpi 21 dpi 
ratio p_ ratio p_ ratio p ratio p _ ratio p ratio p ratio p_ ratio p 
2079, aldehyde dehydrogenase 2 gi 983 00 644918 879 155 0.001 2.72 0.09 
2083 family 118098552 
2751 aldo-keto reductase gi 795 00 367908 7.63 0.99 0.05 248 0.079 
14330324 
2382 alpha-enolase gi 611 00 476175 6.17 487 0.019 2.28 0.91 
46048768 
2481, — annexin A2 (ANXA2) gi 783 00 38901 692 154 0371 N/A 0 244 0.012 4.63 0.003 
4391 45382533 
2224 cartilage matrix protein gi 1090 00 54798.1 6.63 12.82 0.019 259 0 
precursor 71896317 
689 creatine kinase M-type gil211528 = 443 00 43529.1 65 0.87 0.567 5.63 0.001 
4195,  fructose-bisphosphate gil999392 256 00 390228 5.79 0.34 0.009 0.26 0 2.25 0.009 2.23 0.002 
4953 aldolase C 
3273 heat shock protein beta-1 gi 766 00 21715 5.77 542 0.005 2.66 0.955 
(HSPB1) 45384222 
4317 lamin-A gi 651 00 733485 65 288 0 041 0.001 
4310, 45384214 
4321 
2592 nmrA-like family domain- gi 1090 00 32726 =7.77 056 0.001 2.78 0.003 
containing protein 1 71897147 
2483 similar to D4-GDP- gi 552 00 22928.6 5.08 2.37 0.039 1.61 0.254 
dissociation inhibitor 50728568 
(GDID4 
1932 similar to malate gi 819 00 37400.6 8.83 N/A 0.019 13 031 
dehydrogenase 2, NAD 50758110 
(mitochondrial 
5132 similar to myozenin gi 319 00 29589.7 7.85 2.24 0.058 N/A 
50749396 
1919, pyruvate kinase (PKM2) gi212571 291 00 584339 7.29 1.12 0488 044 0.071 05 0.012 255 0.004 281 0.024 1.28 0.088 
1659 
4635 selenium binding protein 1 gi 731 00 530976 6.17 0.82 0.156 1.86 0.034 
118102241 
2016 succinyl-CoA:3-ketoacid- gi 201 00 565494 8.01 549 0.138 2.09 0 
coenzyme A transferase 1, 60592998 
mitochondrial 
4680 UMP-CMP kinase gi 947 00 223863 6.75 1.07 0.56 242 0.002 
71896025 
4786 vimentin (VIM) gil212868 1210 00 53166.6 5.09 478 0 14 0.394 


a) Match ID represents unique numbers assigned to each spot in the matched standard of ImageMaster 2D Platinum version 6.0 software 
b) Accession no. is the MASCOT result of MALDI-TOF/TOF searched from the NCBInr database 

c) The criterion for a successfully identified protein is a protein score confidence interval (C.l. %) for PMF and MS/MS data = 95% 

N/A: A indicates the spot was detectable on one of the gels, N indicates the spot was too weak to detect on one of the gels 
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Table 3 Similar abundance changed proteins in kidney from chicken infected with IBV ck/CH/LDL/97I P, and ck/CH/LDL/97I P,,5 


Master Protein description Accession Protein Protein Mw pl Average ratio and p-value 
ID? no.” score score Cl (Da) 
%S 
P.-infected/control P,15-infected/control 
4 dpi 7 dpi 14 dpi 21 dpi 4 dpi 7 dpi 14 dpi 21 dpi 
ratio p ratio p_ ratio p_ ratio p_ ratio p_ ratio p_ ratio p_ ratio p 
2033 carbonyl reductase [NADPH] _ gi 1110 00 30519.9 85 1.66 0.005 
1 71895267 
435 catalase gi 529 00 60279 809 0.64 0.002 0.67 «20 
53127216 
252, chain A, the structure of gi 989 00 68009.9 655 427 0.015 3.6 0.019 1.56 0.038 
264 chicken mitochondrial Pepck 158430534 
in complex with Pep (PCK2) 
2125, chain A, transthyretin gil1633502 672 00 14209 5.1 0.63 0.005 0.99 088 0.83 0.065 064 0 
3724 
658 chain C, crystal structure of — gi 659 00 47400.8 5.32 171 0.039 1.68 0.022 
native chicken fibrinogen 21730885 
with two different bound 
ligands 
1306 class | alcohol gi 286 00 40891.8 7.85 149 0.001 1.54 0.001 
dehydrogenase, beta subunit 45384164 
2087 cytochrome c oxidase gi 359 00 19676.1 8.91 0.58 0.008 0.57 0.026 
subunit 4 isoform 1, 71895513 
mitochondrial 
1274 D-amino-acid oxidase gi 298 00 41209.8 6.7 14 0.15 1.55 0.001 
118098567 
144 ovotransferrin BB type gi 1020 00 79606.1 685 1.3 0.003 1.56 0 
71274075 
374 similar to CDNA sequence gi 1030 00 655764 785 0.64 0.001 060 O 
BC048390 8098312 
417 similar to CDNA sequence gi 734 00 655764 7.85 0.52 0.008 0.81 0.13 
BC048390 8098312 
422 similar to CDNA sequence gi 865 00 655764 7.85 1.35 0.033 1.52 0.021 
BC048391 8098312 
436 similar to cystathionase gi 618 00 445547 686 0.78 0.032 0.66 0.001 
8094764 
2149 similar to dodecenoyl- gi 502 00 345578 93 1.63 0.003 1.19 0.023 
coenzyme A delta isomerase 8098151 
(3,2 trans-enoyl-coenzyme A 
isomerase) 
670 similar to FLJ20699 protein gi 352 100 528208 5.95 0.75 0.006 0.5 0 
118083181 
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Table 3 Similar abundance changed proteins in kidney from chicken infected with IBV ck/CH/LDL/971 P, and ck/CH/LDL/97I P,,5 (Continued) 


1527, similar to guanidinoacetate gil 424 100 174577 659 063 0 0.65 0.001 0.73 0.001 0.86 0.002 
2150 N-methyltransferase 118103242 
1174 receptor-associated protein gil2661436 440 100 40792.2 8.61 1.66 0.032 1.28 011 
511 retinal dehydrogenase 1 gil 534 100 56400.8 7.49 1.26 0.096 1.57. 0.013 
45383031 

1739 similar to methyltransferase gil 614 100 30640.5 6 1.65 0 

24 50750103 
a) Match ID represents unique numbers assigned to each spot in the matched standard of Ettan™ DeCyder 2D v6.5 software 


b) Accession no. is the MASCOT result of MALDI-TOF/TOF searched from the NCBInr database 
c) The criterion for successful identification of a protein is a protein score confidence interval (C.l. %) for PMF and MS/MS data = 95% 
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Table 4 Differential abundance changed proteins in kidney from chicken infected with IBV ck/CH/LDL/97I1 P; and ck/CH/LDL/97I P,,5 


Master Protein description Accession Protein Protein Mw pl! Average ratio and p-value 
ID? no.P score score Cl (Da) 
%* 
P;-infected/control P,45-infected/control 
4 dpi 7 dpi 14 dpi 21 dpi 4 dpi 7 dpi 14 dpi 21 dpi 
ratio p ratio p_ ratio p_ ratio p_ ratio p_ ratio p_ ratio p_ ratio p 
496 aldehyde dehydrogenase 4. gi 1010 100 112936 54 0.56 0.009 0.85 0.031 
family, member A1 118101121 
662 aldehyde dehydrogenase 4. gi 890 100 112936 54 1.03 084 1.59 0.015 
family, member A1 118101121 
452 alpha-enolase gi 353 100 476175 6.17 125 0.14 17 0.003 
46048768 
821, alpha-tropomyosin of gil833618 94 99.999 32962.7 4.67 2.22 0.003 1.16 046 
1088 smooth muscle (TPM1) 
1977 apolipoprotein A-| (APOA1) giJ211159 = 703 00 30673.2 5.58 0.97 0.31 158 0 
1262 chain A, crystal structures of — gi 728 00 36158.6 5.61 2.04 0.004 11 0.75 
chicken annexin V in 62738641 
complex with Ca?* (ANXAS) 
983, chain C, crystal structure of gil8569623 921 00 47485.9 54 169 0.002 1.24 0.034 1.07 0.35 076 0.006 
784 native chicken fibrinogen 
1926 low molecular weight gi 304 00 18640.2 6.81 6.34 0.045 074 0.16 
phosphotyrosine 86129490 
proteinphosphatase (TCP1) 
1847 manganese-containing gi 258 00 25158.7 86 1.52 0.015 1.12 0.022 
superoxide dismutase 12034955 
precursor (MNSOD) 
1970 nucleoside diphosphate gil2827446 445 00 17542, 7.11 468 0.003 148 081 
kinase 
293 phosphoenolpyruvate gil212538 492 00 702243 6.1 1.15 0.047 1.78 0 
carboxykinase (EC 4.1.1.32) 
665 similar to aflatoxin aldehyde — gi 286 00 369434 6.76 2.61 0.004 1.28 0.059 
reductase 118101125 
789 similar to betaine gi 1010 00 45552.1 7.56 251 141 
homocysteine 50755288 
methyltransferase 
582 glutamate dehydrogenase 1, gi118534 339 00 560754 848 1.74 0.051 1.06 0.17 
mitochondrial 
1461, similar to methyltransferase gi 556 00 30640.5 6 0.72 0.002 1.25 0.074 1.09 0.066 042 0 
2026 24 (MET24) 50750103 
987 sulfotransferase gi 265 00 363335 5.89 453 0.01 0.98 0.85 
45384226 
1650 triosephosphate isomerase gil212774 = 1290 00 26831.9 6.71 1.51 0.009 1.07 043 
(EC 5.3.1.1) 


a) Match ID represents unique numbers assigned to each spot in the matched standard of Ettan™ DeCyder 2D v6.5 software 


b) Accession no. is the MASCOT result of MALDI-TOF/TOF searched from the NCBInr database 


c) The criterion for successful identification of a protein is a protein score confidence interval (C.l. %) for PMF and MS/MS data = 95% 
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Additional file 4. Figures 2 and 3 show representative 
gels indicating protein spots identified in tracheal and 
kidney tissues respectively. As shown in Tables and Fig- 
ures, several proteins were identified in multiple spots 
with the same molecular weight but different isoelectric 
points, such as three spots contained lamin-A (spot # 
4317, 4310, 4321 in Figure 2). 

According to annotations from the UniProt Knowledge- 
base (UniProtKB) and the Gene Ontology databases, most 
of the identified proteins were involved in cytoskeleton 
organization, anti-oxidative stress, the stress response, 
acute phase response, energy metabolism, macromolecular 
biosynthesis, signal transduction and ion transport (sum- 
marized in Additional file 5). Among them, the abundance 
of annexin A2, annexin A5, pyruvate kinase (PKM2), 
alpha-enolase, mitochondrial phosphoenolpyruvate car- 
boxykinase (PCK2), triosephosphate isomerase, heat shock 
protein beta-1, manganese-containing superoxide dismu- 
tase (MnSOD), vimentin, lamin-A, cartilage matrix pro- 
tein, alpha-tropomyosin, nucleoside diphosphate kinase, 
sulfotransferase, and low molecular weight phosphotyro- 
sine proteinphosphatase, were induced to be differentially 
patterns changed in chickens infected with the P; and P1415 
strains, suggested that infection with the Ps and P15 
strains produces different host response. Furthermore, this 
difference in the pattern of change was induced predomi- 
nantly in the early stages of the infection cycle, which sug- 
gests that critical events early in infection are likely to be 
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of key importance in determining the fate of the host. In 
addition, fibrinogen B and y chains showed a similar pat- 
tern of change, with increase in the Ps- and Pj15-infected 
groups. 


Gene ontology annotations of differentially changed 
proteins in abundance 

In order to generate an overview of the subcellular loca- 
tion and biological processes of this abundance changed 
proteins in trachea and kidney, Ps-infected and P,,5- 
infected group, and different dpi, categorization of these 
proteins was performed on the basis of Gene Ontology 
(GO) annotations. As shown in Figure 4A, cellular com- 
ponent ontology revealed that the majority of the identi- 
fied proteins were associated with intracellular 
(GO:0005622) and cytoplasm (GO:0005737) both in tra- 
chea and kidney, P;-infected and P,,5-infected group, 
and different dpi. As shown in Figure 4B, biological pro- 
cess ontology revealed that the majority of the identified 
proteins were associated with metabolic processes 
(GO:0008152), response to stimulus (GO:0050896), and 
regulation of biological process (GO:0050789) both in 
trachea and kidney, P5-infected and Pj,;5-infected group, 
and different dpi. 


Analysis of mRNA level by real-time RT-PCR 
The alterations in the mRNA level of 12 selected pro- 
teins in the trachea or kidney from the control, Ps- 
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Figure 4 Gene Ontology annotation analysis of differentially expressed proteins according to their cellular component (A) and 
biological process (B). This classification was produced on the basis of an analysis using the GOSlimViewer tool at the Agbase database http:// 
www.agbase.msstate.edu/. The results were compared between trachea and kidney, P5-infected and P115-infected group, and different day's 
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infected, and P,,5-infected groups were analyzed at 4, 7, 
14, and 21 dpi. For genes from tracheal tissue, as shown 
in Figure 5A, at 4 dpi, the mRNA level was up-regulated 
in the P;-infected group, when compared with the Pj15- 
infected group. By 7 and 14 dpi, the majority of the 
mRNA exhibited similar down-regulation in both the 
P;-infected group and the P;,5-infected group. By 21 
dpi, their change patterns were different: the mRNA 
level of these genes was up-regulated in the P;-infected 
group, but for the P,,5-infected group, they remained 
low at this time point. For genes from kidney tissue, as 
shown in Figure 5B, at 4 dpi, infection with the P,,5 
strain resulted in a stronger up-regulation in the mRNA 
level than infection with the Ps strain. However, by 7 
dpi, infection with P; induced a stronger down-regula- 
tion than P,,5 infection. At 14 dpi, both virus strains 
resulted in a similar response, with the majority of 
mRNAs showing a down-regulation. In contrast, at 21 
dpi, the majority of mRNAs were still down-regulated in 
the P,,5-infected group, but Ps infection induced up- 
regulation of the mRNA level of several genes. 

In comparison with the results obtained using proteo- 
mics (summarized in Additional file 6 and Additional 
file 7), the trends in the mRNA levels of these genes 
were not completely consistent with the change patterns 
of their corresponding proteins in 2-DE or 2-DIGE gels. 
The disparity between the levels of mRNA and their 
corresponding proteins may occur because posttran- 
scriptional and posttranslational modifications, as well 
as differential mRNA and protein degradation rates, 
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may also contribute to these discrepancies. Several 
papers have addressed this question performing parallel 
proteomic/gene expression studies [23-26]. 


Protein validation by western blot analysis 

To confirm the dynamic alterations of protein abun- 
dance during infection with IBV ck/CH/LDL/97I P; and 
ck/CH/LDL/97I1 Piis, three proteins including annexin 
A2 (ANXAQ2), annexin A5 (ANXAS5), and heat shock 
protein beta-1 (HSPB1) were selected for western blot 
analysis. As shown in Figure 6A and 6B, the abundance 
of HSPB1 was increased in trachea tissues from the P;- 
infected group at 4, 7, and 21 dpi, relative to the control 
group. It was also increased in the Pj;5-infected group, 
relative to the control group. The abundance of Annexin 
A2 was increased in tracheal tissues from the P,,5- 
infected group at 4 and 7 dpi, but decreased in tracheal 
tissues from the Ps-infected group at 7 dpi, relative to 
the control group. In addition, the abundance of 
annexin A5 was increased in Ps-infected kidney tissues 
at 7 dpi, relative to the control group. These data were 
in agreement with the results obtained from the 2-DE 
and 2-DIGE analysis. 


Discussion 

Viruses are obligatory intracellular pathogens that rely 
on the host cell for essentially all steps of their life 
cycle. Although coronaviruses use host proteins as part 
of their replication strategies, it has also become clear 
that the immune, metabolic, stress, cell cycling and 
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different genes, refer to Tables 1, 2, 3, and 4. 


Figure 5 Transcript analysis of the 12 proteins differentially expressed in trachea (A) and kidney (B) of chickens infected with the IBV 
ck/CH/LDL/97I Ps and ck/CH/LDL/97I P115 strains by real-time RT-PCR. Total RNA extracts were prepared from the trachea and kidney of 
chickens in the control, Ps-infected, and P,,5-infected group's at all four time points. Data represent means of three biological replicates per 
group. Error bars indicate standard error. Samples were normalized with the expression of the 18S ribosomal RNA gene. For symbols indicating 
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Figure 6 (A) Western blot analysis of heat shock protein beta-1 (HSPB1) and annexin A2 (ANXA2) in trachea tissue, annexin A5 
(ANXAS5) in kidney tissue of chickens infected with the IBV ck/CH/LDL/97I P; and ck/CH/LDL/97I P,,5 strain. (B) Protein band density was 
analyzed with the BandScan software. Beta-Actin was used as the internal control. Mean values + SE were calculated from three independent 


other pathways are activated by infection. Determina- 
tion, using genomics and proteomics, of the extent to 
which virus-host interaction is coronavirus-specific and 
organ-specific, will be of importance [14,27]. Our pre- 
vious study on the global profile of host protein altera- 
tions in response to IBV infection was focused on an in 
ovo infection model system [21]. In this study, we 
applied a comparative proteomics technical platform for 
the first time to explore the abundance changed protein 
in trachea and kidney tissues from chicken infected in 
vivo with the highly virulent ck/CH/LDL/97I Ps strain 
of IBV and the embryo-passaged, attenuated ck/CH/ 
LDL/97I P,,5 strain of IBV. Some identified proteins are 
likely to be important in the host response to virus 
infection, including cytoskeletal proteins, stress response 
proteins, and anti-oxidative proteins. Interestingly, the 
abundance of these proteins showed different change 


patterns with IBV strains differ in virulence, suggesting 
that some of these differences might be responsible for 
virulence, and consistent with our previous study which 
demonstrated the differences in pathology and virulence 
for these two different viruses [22,28]. These results pro- 
vide an overview of the proteome profile of the host in 
response to different virulent IBV infection in vivo. 

The cytoskeletal network is a cellular scaffold system 
whose functions include maintenance of cellular shape, 
enabling of cellular migration, cell division, intracellular 
transport, signaling, and membrane organization. Some 
host cytoskeletal proteins have been reported as differ- 
entially altered by virus infection in quantitative proteo- 
mic studies [16-18]. Recently, study has revealed that 
several different cellular proteins involved in cell mor- 
phology and the cytoskeleton changed in abundance in 
IBV infected cells [19,20]. IBV infection resulted in a 
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number of changes to the nucleolus both in terms of 
gross morphology and protein content [29]. Our results 
presented a tendency that the abundance of several 
cytoskeletal proteins showed increased with the degree 
of virulence getting strong, which is possibly due to the 
collapse and dispersal of the cytoskeleton in the IBV- 
infected cells, as demonstrated by other viruses [30]. 
Vimentin is a major component of type III intermediate 
filaments that has been reported to be redistributed in 
cells around sites of virus replication and assembly dur- 
ing virus infection. For instance, infection with African 
swine fever virus (ASFV) can lead to the rearrangement 
of vimentin into a cage surrounding a virus factory, 
which may prevent movement of viral components into 
the cytoplasm and concentrate late structural proteins at 
sites of virus assembly [31]. Specifically, vimentin was 
observed to be increased in abundance in the cytoplas- 
mic proteome of IBV-infected cells [19,20]. In our 
study, the abundance of vimentin was increased more 
strongly in the highly virulent Ps strain infected group, 
compared with the attenuated Pj,5 strain infected group 
at 21 dpi. Given that vimentin can confer rigidity to 
domains of the cytoplasm, the vimentin cage may pro- 
vide a physical scaffold to facilitate the construction of 
the virus factory. The stronger induction of vimentin 
could be a result of the large replication capacity of the 
highly virulent IBV strain. Tropomyosin belongs to the 
family of actin-binding proteins that serve important 
functions in microfilament stabilization, regulation of 
microfilament branching, actin polymerization, and 
intracellular transport. The abundance of alpha-tropo- 
myosin was found increased in IBV-infected cell [19,20]. 
In our study, the abundance of alpha-tropomyosin was 
increased much more strongly with highly virulent IBV 
ck/CH/LDL/97I Ps; infection than with the attenuated 
P45 strain, it is likely that it help to orchestrate virus 
assembly, release and efficient cell-to-cell spread, also 
may due to the different virulence between IBV ck/CH/ 
LDL/97I Ps and P,,5 strain. In addition, in our previous 
in ovo infection model [21], the abundance of some 
cytoskeletal including alpha-tropomyosin was decreased 
in the IBV-infected embryonic tissues, this difference 
need to be further investigated. 

In both our previous [21] and current study, the abun- 
dance of annexin A2 and annexin A5 were found to be 
changed upon IBV infection. Annexins are a family of 
conserved proteins characterized by their ability to bind 
and order charged phospholipids in membranes, often 
in response to elevated intracellular calcium. These 
family members are involved specifically in a diverse 
range of cellular functions both inside the cell and 
extracellularly [32]. Annexin A2 (ANXA2) can associate 
with actin filaments and mediates membrane trafficking 
and membrane-cytoskeletal interactions. It has been 
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identified as an important host factor for several viruses 
and at different stages of their life cycle. ANXA2 is 
incorporated into cytomegalovirus and influenza virus 
particles, promotes the entry of virus, and plays a role in 
Human immunodeficiency virus -1 (HIV-1) assembly, 
Bluetongue virus (BTV) release, and hepatitis C virus 
(HCV) replication [33-40]. Annexin A5 (ANXAS5) is 
involved in various intra- and extracellular processes 
including signal transduction, anti-inflammatory pro- 
cesses, membrane trafficking, and ion channel activity 
[41], it also acts to regulate blood coagulation, binding 
to and shielding exposed phospholipids and masking 
their pro-thrombotic properties [42]. The expression of 
ANXA2 and ANXAS5 has been found to be altered dur- 
ing other IBV infection system [19,20]. In this study, the 
abundance of ANXA2 showed weaker increase at early 
stages of infection with the highly virulent IBV ck/CH/ 
LDL/97I Ps strain than with the attenuated P,,5 strain. 
For ANXAS, the level of expression induced at 7 dpi 
with the highly virulent P,; strain of IBV was greater 
than that induced by the attenuated P,,;5 strain. These 
trends were also validated by western blot analysis. 
Whether this difference between P,; and Pj, infectioned 
groups is related with the virulence of IBV or various 
responses of the host cell deserves further study. 

MnSOD is the primary antioxidant enzyme in the 
mitochondria that catalyzes the conversion of superox- 
ide molecules to hydrogen peroxide and molecular oxy- 
gen and therefore forms one of the cell’s major defense 
mechanisms against oxidative stress [43]. Studies have 
reported that viral infections cause oxidative stress, 
which is associated with the activation of phagocytes 
and an increase in the release of reactive oxygen species 
(ROS) that play a positive modulatory role in immune 
activation, the inflammatory response, eradication of 
viral infections and immunity-induced cellular injury 
[44]. The expression of MnSOD was altered during 
infection with some viruses, such as porcine reproduc- 
tive and respiratory syndrome virus, influenza A virus, 
and SARS-CoV [45-47]. The abundance of the anti-oxi- 
dative protein MnSOD was increased much significantly 
in highly virulent IBV ck/CH/LDL/97I Ps infected group 
than the attenuated P,,5 strain infected group. The 
much amplitude increased production of MnSOD dur- 
ing highly virulent strain infection can affect the host 
cell pro-/anti-oxidant balance, which probably results in 
more significant immune activation and a stress 
response induced by reactive oxygen species. It is 
important in the eradication of viral infections and 
immune-induced cellular injury. 

When comparing the protein profiles of groups 
infected with different strains, we found that the abun- 
dance of HSPB1, a well-known heat-shock protein, in 
highly virulent IBV ck/CH/LDL/97I Ps strain infected 
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group was significantly higher than the embryo-pas- 
saged, attenuated P,,5 strain at the early stage of infec- 
tion. Western blot and real-time RT-PCR further 
confirmed this alteration. HSPB1 has different cytopro- 
tective roles, including acting as a molecular chaperone, 
maintaining the normal function of cells through inter- 
action with and stabilization of the cytoskeleton, regula- 
tion of translational initiation, modulation of 
inflammation, inhibition of apoptosis, stimulation of 
innate and adaptive immune responses, and responding 
to a wide variety of stressful stimuli [48-50]. It is prob- 
ably involved in all phases of the viral life cycle, includ- 
ing cell entry, virion disassembly, viral genome 
transcription, replication and morphogenesis. The abun- 
dance of HSPB1 was observed to be increased in IBV- 
infected Vero cells using stable isotope labeling with 
amino acids in cell culture (SILAC) [20]. IBV infection 
induces cell cycle arrest at both S and G2/M phases and 
caspase-dependent apoptosis at late stages of the viral 
infection cycle [51-53]. Results presented in our study 
were probably due to the highly virulent virus triggering 
a markedly more robust inflammation and stress 
response. It also suggests that many important, and 
probably different events in IBV pathogenesis and 
immunology, such as the stress response, inflammation, 
and apoptosis, occur early in infection, and that these 
events may contribute to development of an appropriate 
immune response and the outcome of viral infection. In 
our previous study [21], the abundance of HSPB1 was 
found to be decreased after IBV infection in ovo. This 
difference in different infection model need to be further 
examined. 

In this study, infection with both the highly virulent 
strain and the embryo-passaged, attenuated strain of 
IBV was accompanied by elevation of proteins related to 
energy metabolism. Viruses are obligate parasites that 
are completely dependent on their host’s cellular meta- 
bolism for reproduction. Viral infection has been shown 
to modulate the levels of numerous host metabolic com- 
ponents in pathways such as glycolysis, the tricarboxylic 
acid (TCA) cycle, pentose phosphate pathway, and 
macromolecular biosynthesis in order to redirect valu- 
able resources to their own mass production. Studies 
have demonstrated that the rate of glycolysis in cells 
infected by Rous sarcoma virus, feline leukemia virus, 
and poliomyelitis virus was increased by as much as 
370% [54]. Glycolytic enzymes are involved not only in 
carbohydrate metabolism and increased ATP produc- 
tion, also take part in the induction of anti-oxidative 
stress in host cells and contribute to transcription of 
RNA virus genomes [55]. The abundance of some pro- 
teins involved in glycolysis, such as aldose reductase, 
pyruvate kinase, alpha-enolase, and triosephosphate iso- 
merase was found increased in IBV infected cell [19,20]. 
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Our study also showed that the abundance of some gly- 
colytic enzymes, including phosphoenolpyruvate carbox- 
ykinase (PCK1), mitochondrial phosphoenolpyruvate 
carboxykinase (PCK2), pyruvate kinase, and alpha-eno- 
lase, were increased significantly in the early stage of 
infection with both the highly virulent and the attenu- 
ated IBV strain. These findings support the view that 
IBV infection probably triggers activation of host energy 
metabolic components via systemic or global mechan- 
isms, to keep up with the energy demands of its own 
replication. 

The post-translational modification (PTM) plays criti- 
cal roles in cellular regulation. It has been reported that 
PTM and expression of highly related gene sequences 
can induce multiple protein spots in animals [56]. An 
advantage of 2-DE is the monitoring of multiple forms 
of a protein species, which offers the opportunity to 
investigate the effects of the virus infection on protein 
modification. Our study revealed that several proteins 
were expressed by multiple spots in our study. Although 
we do not know at present how these multiple spots 
were generated, some of them were more likely due to 
post-translational modifications (PTMs), or highly 
related gene sequences coding isoforms [57]. 


Conclusion 

In summary, we investigated the proteome profiles of 
tracheal and kidney tissues from chicken infected with 
highly virulent IBV ck/CH/LDL/97I Ps; and embryo-pas- 
saged, attenuated IBV ck/CH/LDL/97I Pi15. Some pro- 
teins involved in cytoskeleton organization, stress 
response, and anti-oxidative stress, showed different 
change in abundance with IBV strains differing in viru- 
lence. While the roles of some identified alterations 
could be related to host antiviral response or pathogenic 
mechanisms, functional significance of other altered 
proteins remains unclear and needs further investiga- 
tion. However, our findings provide proteome-related 
information on a large scale that should be useful in 
increasing our understanding of the pathogenic and 
immune mechanisms of IBV infection. In addition, they 
will also provide reference for similar research into 
other coronaviruses. 


Materials and methods 

Experimental animals and infecting virus strains 
One-day-old White Leghorn specific-pathogen-free 
(SPF) chickens were obtained from the Laboratory Ani- 
mal Center, Harbin Veterinary Research Institute, the 
Chinese Academy of Agricultural Sciences, China. The 
chickens were maintained in isolators with negative 
pressure, and food and water were provided ad libitum. 
This study was approved by the Animal Welfare Com- 
mittee of Heilongjiang Province, China. 
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The highly virulent ck/CH/LDL/97I Ps strain of IBV 
and attenuated IBV strain ck/CH/LDL/97I P,,5 were 
obtained by passaging the IBV ck/CH/LDL/97I strain, 
which was isolated and identified by our laboratory as 
described previously [22]. The IBV strains were propa- 
gated once in 9- to 11-day-old embryonated SPF 
chicken eggs for preparation of seed stock. The presence 
of viral particles in the allantoic fluid of inoculated eggs 
was confirmed with a negative contrast electronic 
microscope (JEM-1200, EX) as described previously [22]. 


Experimental design 

Sixty-six one-day-old SPF White Leghorn chickens were 
divided randomly into three groups: the Ps-infected 
group, P,,;5-infected group and control group. Twenty- 
two chickens were allocated to each group and housed 
in different isolators. Chickens in the P;-infected and 
Pi,5-infected groups were inoculated with the IBV ck/ 
CH/LDL/97I Ps; strain and ck/CH/LDL/97I Pj15 strain, 
respectively, by oculo-nasal application at 11 days of age 
with a dose of log10° EIDso per chicken. The chickens 
in the control group were mock-inoculated in parallel 
with sterile allantoic fluid. 

Blood samples from 10 chickens in each group were 
collected at 4, 7, 14, and 21 days post inoculation. The 
sera were collected for enzyme-linked immunosorbent 
assay (ELISA) testing. Three chickens were selected ran- 
domly from each group and killed humanely at 4, 7, 14, 
and 21 days post inoculation, respectively. Trachea and 
kidney tissues were separated rapidly and washed with 
ice-cold PBS buffer, snap-frozen in liquid nitrogen, and 
kept subsequently at -80°C for 2-DE or 2-DIGE, real- 
time RT-PCR, and western blot analysis. 


Serum antibody detection 

Serum samples were assayed using a commercial IBV 
antibody test kit (IDEXX Corporation, Westbrook, 
Maine, USA) according to the manufacturer’s instruc- 
tions. Each sample was tested in triplicate. Serum-to- 
positive ratios (S/P ratios) were calculated as described 
previously [2,58]. Individual serum titers were calculated 
from these S/P ratios, evaluated as positive or negative, 
and expressed as OD¢s50 nm Values according to the man- 
ufacturer’s instructions. 


Quantitative analysis of IBV in trachea and kidney by 
real-time RT-PCR 

Viral load of IBV was analyzed using R. M. Jones 
described methods [59]. Tissue samples were ground 
and homogenized, 100 mg tissue homogenates were sus- 
pended in 500 ul phosphate-buffered saline (PBS) con- 
taining 100 yg penicillin and 100 yg streptomycin/ml, 
the suspension were freezed thawing three times, then 
centrifuged at 13,000 x g at 4°C for 5 minutes. RNA 
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was extracted from 200 ul tissue supernatant using TRI- 
zol Reagent (Invitrogen) following the manufacturer’s 
instructions. The real-time RT-PCR assay used the fol- 
lowing primer and probe sequences: IBVF forward pri- 
mer CTA TCG CCA GGG AAA TGT C, IBVR reverse 
primer GCG TCC TAG TGC TGT ACC C, IBV Taq- 
Man® probe FAM-CCTGGAAACGAACGGTA- 
GACCCT-TAMRA [59]. One-step real-time RT-PCR 
reactions were performed using One Step PrimeScript® 
RT-PCR kit (TaKaRa Biotech Co. Ltd., Dalian) on Light- 
Cycler® 480 real-time PCR system (Roche) according to 
the following steps: reverse transcription at 42°C for 10 
min, denaturation at 95°C for 10 s. and 40 cycles with 
95°C for 5 s, 55°C for 20 s, 72°C for 10 s, followed by a 
40°C for 10 s cooling step. All of the samples were 
tested in triplicate in each reaction. The data were ana- 
lyzed using LightCycler® 480 Software Version 1.5. 


Preparation of protein samples 

The frozen trachea or kidney tissues were placed in 
liquid nitrogen and ground thoroughly to a very fine 
powder. Samples of 100 mg tissue powder were dis- 
solved in 500 ul lysing solution containing 7 M urea, 2 
M thiourea, 4% CHAPS, 40 mM DTT, 2% IPG buffer 
pH 3-10, 1% Nuclease Mix and 1% Protease Inhibitor 
Mix (GE Healthcare), then incubated for 2 h at 4°C with 
vortexing once every 15 min, and centrifuged at 15,000 
x g for 1 h at 4°C. The supernatant was collected and 
purified with the PlusOne 2D Clean-up kit (GE Health- 
care). The concentration of each protein sample was 
determined with the PlusOne 2D Quant Kit (GE Health- 
care). Protein samples were aliquoted and stored at -80° 
C for 2-DE or 2-DIGE analysis. 


2-DE analysis of trachea protein samples 

Thirty-six samples of tracheal protein from the three 
groups (Ps-infected, P,,5-infected and control) at 4, 7, 
14, and 21 days post-inoculation (dpi) were analyzed by 
2-DE using the methods previously described [21]. First, 
400 ug protein samples were added to rehydration solu- 
tion (7 M urea, 2 M thiourea, 40 mM DTT, 2% CHAPS, 
0.5% pH 3-10 or pH 4-7 IPG buffer, and 0.002% bromo- 
phenol blue) to make the final volume up to 250 ul, fol- 
lowing Isoelectric focusing (IEF), the IPG strips were 
equilibrated and the second dimension separation was 
performed on 12.5% SDS-polyacrylamide gels on the 
SE600 Ruby system (GE Healthcare). Then, the gels 
were stained with PlusOne Coomassie Blue R-350 (GE 
Healthcare) and scanned with an ImageScanner III (GE 
Healthcare). Quantification analyses were performed 
with Image Master 2D Platinum software v6.0 (GE 
Healthcare). For image analysis, three independent gels 
from the Ps-infected and P,,;5-infected groups were 
compared with those from the corresponding control 
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group at 4, 7, 14, and 21 dpi respectively. The normal- 
ized volume values (vol %) of matched protein spots 
were subjected to Student’s f test using the SPSS statisti- 
cal software package version 16.0. The criterion used to 
define differential expression of spots was that the ratio 
of the vol % in the P;-infected group or the P,5- 
infected group vs. the control group was more than 1.5 
(p < 0.05) or less than 0.67 (p < 0.05). Differentially 
changed protein spots were excised manually from the 
gels and subjected to MS analysis. 


2-DIGE analysis of kidney protein samples 

The kidney protein samples from each group at 4, 7, 14, 
and 21 dpi were used for 2-DIGE analysis respectively. 
The pH of the protein samples was adjusted to 8.5 and 
the protein concentration was adjusted to 5 ug/ul. Equal 
amounts of protein from each sample were pooled 
together as the internal standard. The proteins were 
minimally labeled according to the manufacturer’s 
instructions (CyDye DIGE fluor minimal labeling kit, 
GE Healthcare). The Cy2 was used to label the pooled 
internal standard, and Cy3 and Cy5 were used randomly 
to label samples from the control group, Ps-infected 
group and Pj;5-infected group. To minimize system 
error and inherent biological variation, the sample mul- 
tiplexing was also randomized (see Additional file 8) to 
produce unbiased results. Following the labeling reac- 
tion, 50 ug of each Cy2, Cy3, and Cy5 labeled sample 
was mixed, the pooled sample of each gel was adjusted 
to 450 ul with rehydration buffer (7 M urea, 2 M 
thiourea, 2%CHAPS, 2%v/v IPG buffer,130 mM DTT), 
and loaded subsequently onto 24 cm, linear pH 3-10 
IPG strips (GE Healthcare). Isoelectric focusing (IEF) 
was performed on an Ettan™ IPGphor 3 (GE Health- 
care) using the following protocol: 30 V for 12 h, 200 V 
for 2 h, 500 V for 2 h, 1000 V for 2 h, 8000 V for 3 h, 
and 8000 V for 65000 Vh. Then, the IPG strips were 
equilibrated and the second dimension separation was 
conducted on 12.5% SDS-polyacrylamide gels using an 
Ettan™ DALT six system (GE Healthcare). The CyDye- 
labeled gels were scanned using a Typhoon 9400 (GE 
Healthcare). Image analysis was performed using 
Ettan™ DeCyder Software version v6.5 (GE Healthcare) 
as described in the user’s manual. The statistical analysis 
of changes in protein abundance in different gels was 
performed on the basis of the spot volumes. Protein 
spots with statistically significant results for the Stu- 
dent’s t-test (p < 0.05) and Average Ratio more than1.5 
or less than 0.67 were considered to be differentially 
changed protein spots. 

The preparative gels were made with 1200 ug of unla- 
beled proteins (400 ug from each group). These mixed 
proteins were loaded and separated under the same con- 
ditions as described above. The preparative gels were 
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stained with Coomassie Blue R-350. Each spot of inter- 
est, defined on the basis of CyDye images, was matched 
with a Coomassie Blue R-350 image, then excised 
manually from the gel and subjected to MS analysis. 


MALDI-TOF/TOF MS and database search 

The gel samples were placed in a tube and washed twice 
with 500 ul and 250 ul ddH2O for 15 min. For trypsin 
digestion, the gel samples were washed twice with 50 
mM w/v NH,HCO3 and covered with 10 mg/ml Porcine 
Trypsin solution (Promega, Madison, WI, USA) in 50 
mM w/v NH,HCO3. After incubation overnight at 37°C, 
the supernatant was removed into a second tube and 40 
ul 50 mM w/v NH4HCOs3 was added. Gel samples were 
washed with 40 ul of 50 mM w/v NH4HCOs, the super- 
natant was collected, and both collected supernatants 
were combined. The gel was washed with 70% v/v ACN 
and dried in a Speed Vac (Vacuum Concentrator, Bach- 
hofer). The peptide mixtures were desalted using ZipTip 
C-18 RP tips (Millipore, Billerica, MA, USA) which were 
wetted with 100% ACN and equilibrated with 0.1% TFA. 
Peptide samples, which were redissolved in 10 ml 0.5% 
TFA, were eluted with 50% ACN/0.1% TFA and dried 
in a Speed Vac (Vacuum Concentrator). 

The purified peptides were spotted on a MALDI plate 
and covered with 0.7 ul of 2 mg/ml 3, 5-Dimethoxy-4- 
hydroxycinnamic acid matrix (Sigma) with 10 mM 
NH,H2PO, in 60% ACN. All samples were analyzed by 
MALDI-TOF/TOF MS with a 4700 Proteomics Analyzer 
(Applied Biosystems, Foster City, CA). Monoisotopic 
peak masses were acquired in a mass range of 800-4000 
Da, with a signal-to-noise ratio (S/N) of 200. Five of the 
most intense ion signals, excluding common trypsin 
autolysis peaks and matrix ion signals, were selected as 
precursors for MS/MS acquisition. The peptide mass 
fingerprint (PMF) combined MS/MS data were sub- 
mitted to MASCOT version 3.0 (Matrix Science) for 
identification according to the NCBInr database (release 
16/01/2010, 10343571 sequences, 3528215794 residues). 
The search parameters were set as follows: Gallus, tryp- 
sin cleavage (one missed cleavage allowed), carbamido- 
methylation of cysteine as fixed modification, oxidation 
of methionine as variable modification, peptide mass tol- 
erance set at 100 ppm, fragment tolerance set at 0.8 Da. 
The criterion for successful identification of a protein 
was the protein score confidence interval (C.I. %) = 95%. 


Gene ontology (GO) annotation of differentially expressed 
proteins 

Spot identities were submitted to GORetriever http:// 
www.agbase.msstate.edu/ to obtain the GO annotations. 
If no annotation was returned, GOanna was used to 
retrieve GO annotations assigned on the basis of 
sequence similarities. The resulting annotations were 
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summarized according to the GOA and whole proteome 
GOSlim set using GOSlimViewer [60]. 


Analysis of mRNA levels by real-time RT-PCR 

Specific primers were designed according to the corre- 
sponding gene sequences of the MS-identified proteins 
using Beacon Designer software 7.5 (Primer Biosoft 
International). Information on the primers is listed in 
Additional file 9. Total RNA was extracted using TRIzol 
Reagent (Invitrogen) according to the manufacturer’s 
instructions. Two micrograms of total RNA was reverse 
transcribed with 200 U M-MLV Reverse Transcriptase 
(Invitrogen) and 500 ng Oligo(dT)s as the first strand 
primer in 20 ul reaction solution. Real-time RT-PCR 
was carried out using SYBR® Premix Ex Taq™ II kit 
(Takara) on the LightCycler® 480 real-time PCR system 
(Roche) according to the following steps: 30 s at 95°C, 
40 cycles of denaturation at 95°C for 10 s, and annealing 
and extension at 55°C for 45 s. Each sample was ampli- 
fied in triplicate. Quantitative analysis of the data was 
performed using the 24“ method [61]; samples from 
the control group at 4 dpi were used as the calibrator 
(relative expression = 1), and the 18S ribosomal RNA 
gene was used as an internal reference gene. 


Western blot analysis 

Samples of tracheal and kidney proteins from the P;- 
infected group, P,,5-infected group, and control group 
at all four time points were prepared, and the protein 
concentration was determined as described above. 
Equivalent amounts of total protein were subjected to 
12% SDS-PAGE and then transferred to a nitrocellulose 
membrane. After blocking for 1 h at 37°C, the mem- 
branes were incubated separately with mouse monoclo- 
nal antibody to annexin A5 (sc-32321, Santa Cruz 
Biotechnology, USA), mouse monoclonal antibody to 
HSPB1 (sc-51956, Santa Cruz Biotechnology, USA), and 
goat polyclonal antibody to annexin A2 (sc-1924, Santa 
Cruz Biotechnology, USA) overnight at 4°C. The mem- 
branes were incubated subsequently with IRDye700DX 
conjugated anti-mouse secondary antibody (610-130- 
121, Rockland, Gilbertsville, PA) for 1 h at 37°C, and 
scanned finally on a LI-COR infrared imaging system 
using their Odyssey software (Li-Cor Bioscience, Lin- 
coln, NE). B-actin (sc-47778, Santa Cruz Biotechnology, 
USA) was used as a reference protein to check equal 
loading. Triplicate experiments were performed for each 
sample. Densitometric analysis of protein bands was 
done by using BandScan software (Glyko). 


Additional material 


Additional file 1: Table $1 Serological results post inoculation with 
IBV ck/CH/LDL/971 Ps and ck/CH/LDL/971 P4415. 
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Additional file 2: Figure $1 Summary of changes in protein levels 
over time following infection with IBV ck/CH/LDL/97I Ps and P445. 
The y axis shows the number of differentially expressed protein spots; 
individual spots can be found in Tables 1, 2, 3, and 4. 


Additional file 3: This includes the PMF spectrum and confirmed 
MALDI-TOF-TOF spectrum of differentially expressed protein spots 
in IBV-infected chick tracheal tissues. 


Additional file 4: This includes the PMF spectrum and confirmed 
MALDI-TOF-TOF spectrum of differentially expressed protein spots 
in IBV-infected chick kidney tissues. 


Additional file 5: Table $2 Biological function of differentially 
expressed proteins reported in viral infection. 

Additional file 6: Table $3 Comparison of the fold changes for 
protein abundance observed by 2-DE gel analysis and mRNA level 
obtained by real-time RT-PCR in trachea tissues. 

Additional file 7: Table $4 Comparison of the fold changes for 
protein abundance observed by 2-DIGE gel analysis and mRNA 
expression obtained by real-time RT-PCR in kidney tissues. 


Additional file 8: Table $5 Experiment design of the different 
fluorescent dye labeling. 


Additional file 9: Table S6 The primers of Real-time RT-PCR. 


Abbreviations 

DPI: Days post infection; ACN: Acetonitrile; CHAPS: 3-[(3-cholamidopropyl) 
dimethyl-ammonio] -1-propanesulfonate; DTT: Dithiothreitol; IBV: Infectious 
bronchitis coronavirus; IEF: Isoelectric focusing; IPG: Immobilized pH gradient; 
ALDI-TOF-TOF/MS: Matrix-assisted laser desorption/ionization time-of-flight 
tandem mass spectrometry; PMF: Peptide mass fingerprinting; SDS-PAGE: 
Sodium dodecyl sulfate polyacrylamide gel electrophoresis; SPF: Specific 
pathogen free; TFA: Trifluoroacetic acid; 2-DE: Two-dimensional gel 
electrophoresis; 2-DIGE: Two-dimensional fluorescence difference gel 
electrophoresis. 
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